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In  a  proton  exchange  membrane  fuel  cell  (PEMFC),  flow  field  design  is  an  important  factor  that  influences 
the  distributions  of  current  density  and  water  accumulation.  The  segmented  model  developed  in  prior 
study  is  used  to  investigate  the  effect  of  flow  field  patterns  on  current  density  distribution.  This  model 
predicts  the  distributed  characteristics  of  water  content  in  the  membrane,  relative  humidity  in  the  flow 
channels,  and  water  accumulation  in  the  gas  diffusion  layers  (GDLs). 

Three  single  cells  with  different  flow  field  patterns  are  designed  and  fabricated.  These  three  flow  field 
designs  are  simulated  using  the  segmented  model  and  the  predicted  results  are  compared  and  validated 
by  experimental  data.  This  segmented  model  can  be  used  to  predict  the  effect  of  flow  field  patterns  on 
water  and  current  distributions  before  they  are  machined. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  performance  of  a  PEMFC  is  influenced  by  operating  con¬ 
ditions  such  as  current  and  temperature,  and  cell  designs  [1,2]. 
Among  the  fuel  cell  stack  operating  conditions,  temperature,  flow 
rates,  pressure,  and  relative  humidity  (RH)  of  reactants  are  key 
factors  that  affect  membrane  hydration  level,  and  thus  the  elec¬ 
trical  conductivity  of  the  membrane.  As  to  the  cell  designs,  flow 
field  patterns  on  the  bipolar  plates  play  an  important  role  in  dis¬ 
tributing  reactants  evenly  throughout  the  cell,  and  avoid  flooding 
by  removing  generated  liquid  water. 

A  PEMFC  stack  consists  of  multiple  cells  connected  by  bipolar 
plates.  They  account  for  the  majority  of  stack  volume  and  weight. 
The  bipolar  plates  play  multiple  roles  such  as  current  distribution 
and  collection,  transporting  and  distributing  reactants,  and  remov¬ 
ing  generated  liquid  water.  The  power  density  of  a  PEMFC  can  be 
significantly  influenced  by  the  flow  channel  patterns  on  the  bipolar 
plates.  An  appropriate  flow  field  design  can  increase  the  power  den¬ 
sity  and  improve  water  and  heat  managements.  In  some  designs, 
cooling  channels  are  integrated  on  the  bipolar  plates  [3,4].  Thus, 
the  flow  field  pattern  on  the  bipolar  plates  is  an  important  factor  in 
the  design  of  a  PEMFC  stack. 


*  Corresponding  author.  Tel.:  +886  2  26712711x313860;  fax:  +886  2  26721134. 
E-mail  address:  songchen@umich.edu  (Y.-S.  Chen). 


Common  flow  field  patterns  include  parallel,  serpentine,  inter- 
digitated  as  well  as  their  combinations.  Um  et  al.  [5]  used 
computational  fluid  dynamics  (CFD)  method  to  develop  a  three- 
dimensional  model  to  study  the  performance  of  an  interdigitated 
flow  field  design.  Their  results  show  improved  performance  of  fuel 
cell  due  to  forced  convection  of  gases  through  the  gas  GDL.  Wang  et 
al.  [6]  and  Pasaogullari  and  Wang  [7]  developed  a  two-phase  flow 
model  for  the  cathode  GDL,  which  successfully  described  water 
vapor  distribution  and  liquid  water  accumulation  in  the  GDL  and 
in  the  flow  channel.  Due  to  heavy  computation  load,  CFD  models 
usually  have  difficulties  simulating  a  fuel  cell  with  a  complicated 
flow  field  pattern.  Thus,  the  published  CFD  models  typically  were 
demonstrated  using  a  fuel  cell  with  simple  flow  field  such  as  a 
straight  flow  channel. 

Many  studies  have  compared  the  performance  of  fuel  cells  with 
different  flow  field  designs.  Arico  et  al.  [8  ]  conducted  an  experiment 
to  compare  the  influence  of  the  flow  field  design  on  the  perfor¬ 
mance  of  direct  methanol  fuel  cells  (DMFCs).  Dohle  et  al.  [9]  and 
Bewer  et  al.  [10,11]  also  investigated  the  influence  of  flow  field 
design  on  the  performance  of  PEMFCs.  They  concluded  that  flow 
field  patterns  can  significantly  affect  the  distribution  of  gas  flow 
and  current  density  of  a  PEMFC. 

Kumar  and  Reddy  [12]  developed  a  three-dimensional  single 
phase  isothermal  model  to  describe  the  steady-state  and  tran¬ 
sient  response  of  four  PEMFCs  with  different  flow  field  designs, 
including  serpentine,  parallel,  multi-parallel,  and  interdigitated 
types.  Their  modeling  results  showed  that  transient  responses  of 
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Nomenclature 

a  water  vapor  activity 

F  Faraday’s  constant  (96,485  C  equivalent-1 ) 

I  current  (A) 

N  molar  flow  rate  (mol  s-1) 

P  pressure  (Pa) 

V  voltage  (V) 

Greek 

X  water  content 

Subscripts 
an  anode 

avg  average 

ca  cathode 

cell  cell 

H2  hydrogen 

in  inlet 

N2  nitrogen 

out  outlet 

02  oxygen 

pern  proton  exchange  membrane 
seg  segment 

w  liquid  water 


serpentine  and  parallel  designs  were  faster  but  performances  were 
lower  than  the  other  two  designs.  However,  Kumar’s  study  did  not 
consider  water  accumulation  in  the  GDLs.  Due  to  heavy  computa¬ 
tional  method,  their  model  was  not  able  to  simulate  PEMFCs  with 
different  anode  and  cathode  flow  field  designs. 

Lumped  models  which  assume  uniform  reactions  within  the 
whole  area  of  the  fuel  cells  do  not  consider  the  spatial  distribution 
of  reactants.  Golbert  and  Lewin  [13]  developed  a  dynamic  fuel  cell 
model,  which  accounts  for  spatial  dependencies  of  voltage,  current, 
reactant  flows,  and  temperatures  in  the  flow  channel.  However, 
their  model  can  only  describe  the  transient  behavior  within  a  sin¬ 
gle  straight  flow  channel.  In  practice,  the  fuel  cell  performance  and 
water  distribution  depend  on  flow  field  design,  and  the  water  gen¬ 
eration  and  accumulation  inside  the  fuel  cells  distribute  unevenly 
along  the  flow  channels.  A  model  that  can  describe  the  distributed 
properties  of  a  single  cell,  especially  with  complex  flow  pattern,  is 
needed. 

In  order  to  solve  this  problem,  in  our  previous  work  [14]  we 
developed  a  segmented  model  to  predict  distributions  of  liquid 
water  accumulation,  current  density,  water  content  in  the  mem¬ 
brane,  and  relative  humidity  in  the  flow  channels.  In  this  study,  we 
use  the  segmented  model  to  study  four  single  cells  with  different 
flow  field  patterns.  The  performance  and  distributed  properties  of 
these  single  cells  under  various  working  conditions  are  simulated 
and  discussed. 


used  to  validate  modeling  results.  Based  on  the  modeling  results 
of  these  three  patterns,  this  model  approach  was  found  to  pro¬ 
duce  qualitatively  accurate  predictions.  Subsequently,  we  designed 
a  fourth  flow  pattern.  The  modeling  results  of  these  four  cells  are 
then  analyzed  and  discussed.  This  design  procedure  is  depicted  in 
Fig.  2. 

2.1.  Summary  of  the  segmented  model 

In  this  paper,  we  use  the  segmented  model  developed  in  our  pre¬ 
vious  paper  [14]  to  investigate  the  influence  of  flow  field  patterns 
on  cell  performance.  The  details  of  the  model  have  been  published; 
therefore,  only  the  concept  is  summarized  below. 

To  capture  the  distributed  behavior  inside  a  PEMFC,  the  active 
area  is  divided  into  fifteen  segments  and  each  segment  is  viewed 
as  a  small  lumped  model.  These  segments  are  connected  together 
according  to  the  flow  field  directions  of  the  anode  and  cathode, 
as  shown  in  Fig.  3.  Each  segment  consists  of  six  interacting  sub¬ 
models:  cathode  flow  channel,  anode  flow  channel,  cathode  GDL, 
anode  GDL,  membrane  hydration,  and  segment  voltage.  The  model 
is  implemented  in  MATLAB/SIMULINK.  The  block  diagrams  of  these 
six  sub-models  of  a  segment  are  shown  in  Fig.  4.  Given  the  inflow 
properties  and  (initial  guess  of)  a  segment  current,  each  segment 
model  calculates  its  segment  voltage.  Since  the  current  density  is 
actually  not  uniform  throughout  the  active  area,  the  calculated 
segment  voltages  may  not  be  the  same.  Thus,  the  cell  voltage  is 
determined  iteratively  by  the  process  shown  in  Fig.  5.  This  algo¬ 
rithm  is  explained  briefly  below. 

At  the  beginning  of  an  iteration,  we  guess  an  inlet  gas  pressure 
Pin  for  a  certain  load  current  Ice\\.  Under  desired  stoichiometry  val¬ 
ues  and  relative  humidity  levels,  the  inflow  anode  and  cathode 
flow  properties  for  the  cell  are  determined.  The  initial  guess  for 
each  segment  current  Jsegi,  is  set  to  be  one  fifteenth  of  the  cell 
current,  Ice\\.  The  segment  voltage  Vseg  i  is  determined  by  the  seg¬ 
ment  model.  Because  of  different  locations,  each  segment  may  not 
have  the  same  inflow  gas  properties  such  as  reactant  pressure,  rel¬ 
ative  humidity,  and  reactant  concentration,  resulting  in  different 
segment  voltages— which  in  theory  should  be  identical. 

If  the  difference  between  the  maximum  and  minimum  segment 
voltages  is  not  within  an  acceptable  range,  the  segment  currents 
need  to  be  corrected.  Based  on  the  typical  polarization  curve  of 
a  fuel  cell,  the  segment  with  higher  voltage  should  have  higher 
current,  and  the  one  with  lower  voltage  should  lower  its  current. 
To  increase  iteration  speed,  the  change  in  segment  current  is  set 
to  be  proportional  to  the  voltage  difference  while  keeping  the  cell 
current  constant. 

AJi  =  k[max(Vrseg  j)  —  min(Vrseg  j)]  (1) 


/  .1*1+1  r  in, _ fri _ 

seg.i  seg.i  -r  number  0f  segments  with  Vseg  i  >  Vavg 

when  Vseg>i>VaVg  (2) 


2.  Modeling  procedure 

In  this  work,  we  design  three  single  cells  with  different  cathode 
flow  field  patterns,  as  shown  in  Fig.  1 .  These  three  designs  represent 
popular  straight  channel,  serpentine,  and  perpendicular/counter¬ 
flow  flow  field  arrangements.  Three  models  corresponding  to  these 
three  patterns  are  also  developed  to  study  their  influence  on  cell 
performance  and  water  accumulation.  The  model  calibration  is 
done  based  on  experimental  data  of  the  single  cell  No.  1.  Then  we 
keep  the  calibrated  parameters  and  reconnect  the  segments  within 
the  models  to  simulate  the  flow  field  patterns  of  cells  No.  2  and  No. 

3,  respectively.  The  experimental  results  of  cells  No.  2  and  No.  3  are 


segj  segj  number  of  segments  with  l/segj  <  l/avg 

when  Vseg>j  <  Vavg  (3) 

The  superscript  n  in  these  equations  indicates  the  n-th  step  in 
time.  In  addition,  the  inlet  pressures  of  the  anode  and  cathode  are 
also  adjusted  to  keep  the  outlet  pressure  the  same  as  the  ambient 
pressure.  At  the  end  of  the  iteration  procedure,  averaged  values  in 
each  segment  are  calculated  and  assumed  to  represent  the  variable 
values  at  the  center  of  each  segment.  These  values  are  then  used  to 
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(a)  Cell  No.  1 
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(c)  Cell  No.  3 


Fig.  1.  Flow  field  designs  in  this  study  (Left:  anode;  Right:  cathode.),  (a)  Cell  No.  1 ;  (b)  cell  No.  2;  (c)  cell  No.  3. 


Fig.  2.  Flow  diagram  showing  modeling  procedure. 
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Fig.  3.  Schematic  of  a  single  cell  modeled  as  fifteen  segments. 


calculate  the  distributed  values  of  the  whole  active  area  of  the  fuel 
cell  through  interpolation  and  extrapolation. 

2.2.  Modeling  of  counter-flow  and  co-flow  designs 

In  the  flow  field  patterns  of  cells  No.  2  and  No.  3,  counter¬ 
flow  exists  in  part  of  the  active  area,  as  shown  in  Fig.  6.  Fig.  7 
shows  segments  numbered  successively  according  to  the  anode 
flow  direction.  If  we  follow  the  iteration  procedure  shown  in  Fig.  5, 
it  will  be  a  challenge  in  solving  the  counter-flow  segments.  In  this 
study,  segment  current  and  all  the  properties  of  the  anode  and 
cathode  reactants  are  required  for  individual  segment  to  obtain 
the  segment  voltage.  In  the  model  of  cell  No.  2,  if  we  solve  each 
segment  successively  according  to  the  anode  flow  direction,  then 


there  is  no  problem  to  get  Nan  H2  in  and  Nan  win  for  segments  1,  2, 
and  3.  Flowever,  when  we  start  to  solve  segment  4,  we  do  not  have 
the  values  of  Nca  o2,in>  Nca  N2iin,  and  Ncawin  since  they  are  outputs  of 
segment  5,  which  has  not  been  solved  yet. 

A  few  researchers  have  modeled  the  counter-flow  configuration 
of  PEMFCs  [15-18].  Ge  and  Yi  [15]  developed  a  two  dimensional 
model  to  investigate  the  co-flow  and  counter-flow  modes  at  differ¬ 
ent  operating  conditions.  When  they  dealt  with  the  counter-flow 
mode,  they  guessed  a  water  flow  rate  at  the  cathode  inlet  and  the 
calculated  water  flow  rate  at  the  cathode  outlet  was  compared  with 
the  actual  value.  If  these  two  values  do  not  match,  they  guess  a  new 
value  and  iterate  through  the  process  until  convergence.  Birgers- 
son  and  Vynnycky  [16]  developed  an  isothermal  three-dimensional 
model  to  compare  the  performance  of  four  different  flow  field  pat- 


Fig.  4.  The  connection  of  the  six  sub-models  of  a  segment  in  the  SIMULINK  environment. 
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Fig.  5.  Flow  diagram  of  the  iterative  procedure  to  calculate  cell  voltage. 


terns:  co-flow,  counter-flow,  interdigitated,  and  foam.  The  details 
of  solving  counter-flow  mode  were  not  given  in  their  study.  How¬ 
ever,  their  model  took  one  to  twelve  hours  to  solve  the  counter-flow 
model  problem.  Their  results  suggested  that  under  the  same  oper¬ 
ating  conditions,  interdigitated  and  foam  flow  fields  showed  higher 
performance  than  co-flow  and  counter-flow  did,  but  there  was 
no  significant  performance  difference  between  the  co-flow  and 
counter-flow.  Hwang  et  al.  [17]  developed  a  non-isothermal  model 


to  investigate  an  interdigitated  flow  field  pattern.  They  studied 
the  conditions  where  the  cathode  reactants  flowed  into  the  GDL 
with  the  same  (co-flow)  and  opposite  (counter-flow)  directions  to 
the  anode  reactants.  They  used  commercial  software  to  solve  their 
model  and  the  details  were  not  provided  in  their  paper.  Their  model 
also  took  hours  to  solve  for  each  operating  condition. 

Berger  et  al.  [18]  developed  a  simplified  model  to  study  the  dis¬ 
tributions  of  current  density  and  relative  humidity  along  a  straight 
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Fig.  6.  Schematic  representation  of  counter-flow  and  co-flow  modes  in  cells  No.  2  and  No.  3.  (a)  Cathode  of  cell  No.  2;  (b)  anode  of  cell  No.  2;  (c)  cathode  of  cell  No.  3;  (b) 
anode  of  cell  No.  3. 


channel.  They  discussed  the  influences  of  co-flow  and  counter-flow 
modes  on  the  distributions  of  current  density  and  relative  humidity. 
In  their  model  they  discretized  the  flow  channel  into  1048  points. 
When  they  solved  the  counter-flow  mode,  all  the  initial  values 
of  the  anode  and  cathode  fluxes  at  these  points  were  set  as  the 
inlet  values.  The  cathode  fluxes  were  then  updated  by  using  the 
anode  fluxes  from  the  previous  iteration.  The  anode  fluxes  were 
then  updated  by  using  local  current  and  water  transfer  calculated 
from  previously  computed  cathode  channel  values.  The  updating 
process  is  continued  until  the  anode  and  the  cathode  flux  values 


converge.  The  details  can  be  found  in  Ref.  [18].  This  is  a  useful 
method  for  solving  the  counter-flow  problem  in  a  straight  chan¬ 
nel  design.  In  our  study,  there  are  not  only  counter-flow  modes  but 
also  co-flow  modes  in  our  designs  of  cells  No.  2  and  No.  3 ;  moreover, 
the  flow  field  patterns  are  not  simply  straight  channels.  Therefore, 
we  used  a  modified  version  of  Berger’s  updating  procedure  to  solve 
the  problem  in  our  models. 

As  mentioned  at  the  beginning  of  this  section,  Nca  o2,in*  NcaN2iin, 
and  Nca  win  for  segment  4  are  unknown  when  we  start  to  solve 
segment  4  in  cell  No.  2  during  each  iteration.  However,  we  can 


(a) 


(b) 


Anode  flow 


Cathode  flow 


Fig.  7.  Segment  connection  of  cell  No.  2  (a)  and  No.  3  (b). 
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calculate  Nca  o2,in  and  Nca  N2  in  before  segment  5  and  segment  6  are 
solved.  Flow  rates  of  reacted  oxygen  in  segments  5  and  6  can  be 
determined  by  the  assigned  segment  currents  at  the  beginning  of 
each  iteration. 

M  -I  -  M  -I  _  Iseg’5  -  Iseg’6  (A) 

ca,02,in|seg4  ca,02,in|seg6  4/74/7  ^  ' 

Besides,  nitrogen  does  not  react,  so  nitrogen  flow  rate  in  seg¬ 
ment  4  is  equal  to  that  at  the  inlet  of  segment  6. 

Nca,N2,in|seg?4  =  NCa,N2,in  |  seg6  (5) 

Accordingly,  the  only  unknown  is  Ncawin  when  we  start  to  solve 
segment  4.  In  the  flow  field  design  of  cell  No.  2,  same  situation 
happens  to  segments  5, 10,  and  11. 

Similar  to  Berger’s  iteration  procedure,  in  the  beginning  of  our 
solution  process,  the  Nca  w  in  values  of  segments  4, 5, 10,  and  1 1  are 
initially  set  as  the  inlet  values  of  segments  6  or  12  based  on  the 
flow  field  pattern.  Then  Nca>win  values  of  the  above  four  segments 
are  updated  by  using  the  output  values  of  their  preceding  segments, 
i.e.,  segments  5,  6, 11, 12,  respectively. 

JW  ln+1  at  \n 

Nca,w(in|seg4=  NCa,w,out|seg5  (6) 

The  superscript  n  in  these  equations  indicates  the  n-th  step  in 
time.  The  updating  procedure  continues  until  solutions  converge. 
The  same  solution  procedure  is  applied  to  cell  No.  3  for  those  seg¬ 
ments  with  counter-flow  mode. 

The  modeling  results  will  be  discussed  in  the  next  section.  After 
this  model  is  calibrated  and  validated  by  experimental  data,  it  can 
be  applied  to  other  cells  with  different  flow  field  designs. 

3.  Modeling  results  and  discussions 

3.1.  Modeling  results  of  cells  No.  1, 2,  and  3 

Fig.  8  shows  performance  curves  of  cells  No.  1 , 2,  and  3.  All  three 
cells  operated  at  low  cathode  inlet  RH  of  50%  show  lower  cell  perfor¬ 
mance  than  at  inlet  RH  of  100%.  In  addition,  cell  No.  1  shows  lower 
performance  than  cells  No.  2  and  3  when  the  cathode  inlet  RH  is 
50%.  This  is  because  of  the  low  water  content  in  the  cathode  gas 
flow,  and  in  the  design  of  cell  No.  1,  more  water  seems  to  be  taken 
from  the  membrane,  resulting  in  low  membrane  hydration  and  low 
conductivity,  as  shown  in  Fig.  8(b)  and  (c).  The  experimental  results 
suggest  that  cathode  inlet  RH  has  more  significant  influence  on  cell 
performance  in  cell  No.  1 .  A  possible  reason  is  because  the  counter¬ 
flow  mode  exists  in  both  cells  No.  2  and  3  which  helps  to  improve 
membrane  hydration  and  conductivity. 

For  the  model  developed  to  describe  cell  No.  1,  the  parameters 
were  calibrated  based  on  the  experimental  data.  Fig.  8(a)  shows 
that  simulation  results  agree  well  with  the  experimental  results 
after  the  model  calibration.  In  Fig.  8(b)  and  (c),  the  slight  deviation 
between  modeling  results  and  experimental  results  for  cells  No.  2 
and  No.  3  could  indicate  that  the  model  calibrated  well  for  one  flow 
field  design  does  not  fully  capture  the  cell  with  different  cathode 
flow  field  designs,  even  after  we  reconnect  the  15  segments  prop¬ 
erly.  Another  possible  reason  of  this  slight  mismatch  is  that  the 
MEAs  used  for  the  three  single  cells  are  not  identical.  For  example, 
the  difference  of  Pt  catalyst  loading  between  MEAs  could  result  in 
slightly  different  cell  performance.  In  addition,  different  flow  field 
designs  could  also  change  contact  resistances.  These  problems  can 
be  solved  by  tuning  some  parameters  in  the  cell  voltage  model. 
Overall,  we  judge  the  model  prediction  accuracy  to  be  quite  sat¬ 
isfactory.  We  think  this  segmented  model  can  be  reconnected  to 
predict  the  performance  of  cells  with  different  flow  field  design, 
and  this  is  useful  as  a  tool  to  assess  cell  designs  before  those  cells 
are  actually  machined,  assembled  and  tested. 


Fig.  8.  Comparison  of  experimental  results  and  modeling  results  for  different  cath¬ 
ode  inlet  RH  values,  (a)  Cell  No.  1 ;  (b)  cell  No.  2;  (c)  cell  No.  3. 


Fig.  9  shows  current  density  distribution,  water  content  in  the 
membrane,  and  RH  values  in  the  anode  and  cathode  channels  of 
cell  No.  2  when  the  cathode  inlet  RH  is  100%  and  50%,  respectively. 
When  cathode  inlet  RH  is  1 00%,  current  density  distribution  is  more 
even  than  the  case  when  the  cathode  inlet  RH  is  50%.  When  cathode 
inlet  RH  is  50%,  the  current  density  near  cathode  inlet  is  lower  and 
increases  with  cathode  gas  flow  distance,  as  shown  in  Fig.  9(b).  The 
results  also  indicate  that  the  cathode  inlet  RH  significantly  influ¬ 
ences  current  density  distribution.  The  current  density  distribution 
of  cell  No.  3  shows  similar  trend  as  that  of  cell  No.  2,  i.e.,  the  current 
density  near  the  cathode  inlet  is  the  lowest  and  gradually  increases 
toward  the  outlet,  as  shown  in  Fig.  10(b). 

3.2.  Comparison  between  modeling  and  experimental  results 

Figs.  11-13  compare  the  experimental  data  from  RH  sensors 
[19]  and  modeling  results  of  RH  distribution  in  the  flow  channels 
for  cell  No.  3.  Fig.  1 1(c)  shows  that  at  high  load,  anode  RH  near  the 
downstream  area  is  lower  than  that  at  low  current  load,  as  shown  in 
Fig.  11(b).  A  possible  reason  is  the  higher  dry  hydrogen  flow  rate  at 
high  current  density,  resulting  in  decreased  RH  in  the  down  stream. 
The  experimental  data  supports  the  modeling  results,  as  shown  in 
Fig.  11(c).  In  Fig.  1 2,  when  the  cathode  inlet  RH  is  50%,  the  modeling 
results  also  show  good  agreement  with  experimental  data. 

In  Fig.  13(b),  at  low  load  the  RH  in  the  cathode  channels  also 
shows  similar  results  as  experimental  data.  However,  at  high  cur¬ 
rent  density  condition,  the  model  predicted  RH  in  the  downstream 
area  is  100%,  as  shown  in  Fig.  13(c),  but  test  results  are  around  60% 
as  presented  in  Fig.  13(a).  A  possible  reason  is  that  we  assumed 
that  in  our  model  liquid  water  is  moved  out  of  the  cell  at  the 
same  velocity  as  gas,  so  there  is  some  liquid  water  in  the  down¬ 
stream  area  of  the  cathode  channels.  However,  because  channel 
number  is  reduced  from  twelve  to  six,  gas  velocity  in  the  down¬ 
stream  channels  increases,  resulting  in  the  quick  removal  of  liquid 
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Fig.  9.  Distributions  of  current  density,  water  content  in  the  membrane,  RH  in  the  anode  channels, 
right  figures:  cathode  inlet  RH  =  50%  (cell  No.  2). 


and  RH  in  the  cathode  channels.  Left  figures:  cathode  inlet  RH  =  100%; 


water.  Besides,  the  under-saturated  cathode  gas  also  contributes  to 
lower  RH  in  the  downstream  channels. 

Figs.  14  and  1 5  compare  modeling  results  and  experimental  data 
of  water  distribution  in  cells  No.  2  and  3,  respectively.  In  Fig.  14, 
modeling  results  agree  well  with  experimental  data  for  most  of 
the  operating  conditions  except  at  low  current  loads.  According 
to  the  results  shown  in  Fig.  15,  this  model  successfully  predicts 
water  accumulation  in  the  majority  of  segments  for  most  of  the 
operating  conditions.  The  variation  between  modeling  results  and 
experimental  data  could  be  due  to  a  few  different  reasons.  First, 
the  model  parameters  were  calibrated  based  on  the  experimen¬ 
tal  data  of  cell  No.  1.  The  influence  of  flow  field  patterns  on  the 
parameters  was  not  considered.  Second,  the  experimental  data 


of  water  accumulation  was  quantified  “indirectly”  by  using  the 
least-squares  method,  which  was  described  in  our  previous  study 

[1 9] .  An  improved  experimental  design  to  “directly”  measure  water 
accumulation  in  the  GDL  is  needed  for  accuracy  when  water  accu¬ 
mulation  is  considered  in  the  model. 

3.3.  New  flow  field  design 

According  to  the  above  results,  high  RH  value  of  cathode  inlet 
gas  helps  the  cell  current  to  be  distributed  more  evenly.  Liu  et  al. 

[20]  measured  the  current  density  distribution  of  a  working  fuel 
cell;  they  also  observed  this  phenomenon  in  their  study.  This  phe¬ 
nomenon  can  be  explained  by  the  membrane  hydration  model,  in 
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Fig.  10.  Distributions  of  current  density,  water  content  in  the  membrane,  RH  in  the  anode  channels,  and  RH  in  the  cathode  channels.  Left  figures:  cathode  inlet  RH  =  100%; 
right  figures:  cathode  inlet  RH  =  50%  (cell  No.  3). 


which  water  content  in  the  membrane  is  determined  by  the  RH 
values  in  the  anode  and  cathode  [21  ],  as  shown  below: 


A,pem — 0.043  +  1 7.81  cipem- 39.85Qnprn+36.0fll 


Qpem  — 


pern- 


^pem  ' 


apem  > 


Qan  +  Qca 


0  <  Qpem  <  1 

(10) 

(11) 


At  the  operating  conditions  where  the  cathode  inlet  RH  is  50%, 
the  RH  values  in  the  anode  and  cathode  channels  increase  with  flow 
direction.  If  the  cathode  flow  direction  is  opposite  to  the  anode  flow 
direction,  the  overall  water  content  in  the  membrane  will  be  more 
even  due  to  the  fact  reactants  are  more  humid  at  the  downstream. 
When  a  humidifier  is  not  available,  counter-flow  field  design  is  a 
good  solution. 


In  portable  fuel  cell  systems  or  transportation  applications, 
where  there  is  limited  space  for  humidifiers,  an  appropriate  flow 
field  design  may  help  to  increase  the  system  efficiency.  Under  the 
circumstances  where  humidifiers  cannot  work  efficiently  or  are  not 
available,  generated  liquid  water  within  the  fuel  cell  should  be  used 
effectively  to  hydrate  the  membrane  internally.  In  the  flow  pattern 
design  of  cells  No.  1,  2,  and  3,  both  the  anode  and  cathode  reac¬ 
tants  move  from  the  top  to  the  bottom.  Thus,  we  designed  a  virtual 
cell  No.  4,  in  which  the  cathode  inlet  is  placed  at  the  bottom  of  the 
active  area  and  the  anode  flow  direction  is  kept  the  same  as  the 
other  cells,  as  shown  in  Fig.  16.  The  performance  and  water  accu¬ 
mulation  in  this  design  is  predicted  by  our  segmented  model  and 
compared  with  other  designs. 

Fig.  17  shows  the  modeling  results  of  cell  No.  4.  In  the  working 
condition  where  the  cathode  inlet  RH  is  100%,  the  distributions  of 
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Fig.  11.  RH  in  the  anode  flow  channels  when  cathode  inlet  RH  is  100%  (cell  No.  3)  (a)  experimental  results;  (b)  modeling  results  when  current  density  is  0.1  Acm-2;  (c) 
modeling  results  when  current  density  is  0.7  Acm-2. 


current  density  and  water  content  in  the  membrane  are  similar 
to  those  of  cells  No.  1,  2,  and  3  at  the  same  operating  condition. 
However,  when  the  cathode  inlet  RH  is  50%,  the  current  density 
near  the  anode  inlet  of  cell  No.  4  (4.72%)  is  higher  than  those  of 
cells  No.  1  (3.90%)  [14],  No.  2  (3.59%),  and  No.  3  (3.75%).  The  current 
densities  at  the  bottom  segments  of  cell  No.  4  are  lower  because 
under-saturated  cathode  reactant  is  fed  from  the  bottom  of  the  cell. 
The  maximum  current  density  appears  in  the  middle  part  of  the 
active  area,  as  shown  in  Fig.  17(b),  and  the  distribution  is  more 


uniform  through  the  whole  active  area  when  compared  with  that  of 
cells  No,  1,2,  and  3.  Water  content  in  the  membrane  shows  similar 
trend  as  current  density  distribution,  as  shown  in  Fig.  17(d).  The 
performance  of  these  four  flow  pattern  designs  was  compared  in 
Fig.  18.  There  is  no  significant  difference  when  the  cathode  inlet 
RH  is  100%.  When  the  RH  of  cathode  inlet  is  50%,  however,  cell  No. 
4  shows  higher  performance  than  the  other  cells  do,  especially  at 
high  load.  This  is  mainly  due  to  the  humid  cathode  reactant  near 
the  outlet,  resulting  in  increased  water  content  in  the  membrane. 
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Fig.  12.  RH  in  the  anode  flow  channels  when  cathode  inlet  RH  is  50%  (cell  No.  3)  (a)  experimental  results;  (b)  modeling  results  when  current  density  is  0.1  Acm-2;  (c) 
modeling  results  when  current  density  is  0.6  Acm-2. 
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Fig.  13.  RH  in  the  cathode  flow  channels  when  cathode  inlet  RH  is  50%  (cell  No.  3)  (a)  experimental  results;  (b)  modeling  results  when  current  density  is  0.1  Acm-2;  (c) 
modeling  results  when  current  density  is  0.6  Acm-2. 
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Fig.  14.  Comparison  of  average  liquid  saturations  from  modeling  results  and  experimental  data  (cell  No.  2)  (a)  distribution  of  average  liquid  saturation  in  the  cathode  GDL 
under  rib;  (b)  distribution  of  average  liquid  saturation  in  the  cathode  GDL  under  channel. 
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Fig.  15.  Comparison  of  average  liquid  saturations  from  modeling  results  and  experimental  data  (cell  No.  3)  (a)  distribution  of  average  liquid  saturation  in  the  cathode  GDL 
under  rib;  (b)  distribution  of  average  liquid  saturation  in  the  cathode  GDL  under  channel. 
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Fig.  16.  Flow  field  pattern  of  cell  No.  4.  (a)  The  cathode  inlet  is  placed  at  the  bottom  of  the  active  area,  (b)  The  anode  inlet  is  on  the  upper  right  corner. 
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Fig.  17.  Distributions  of  current  density,  water  content  in  the  membrane,  RH  in  the  anode  channels,  and  RH  in  the  cathode  channels.  Left  figures:  cathode  inlet  RH  =  100%; 
right  figures:  cathode  inlet  RH  =  50%  (cell  No.  4). 


Fig.  18.  Comparison  of  performances  of  four  single  cells. 


4.  Conclusions 

In  this  work  current  density  distributions  within  four  single 
cells  with  different  flow  field  patterns  are  predicted  by  using  a  seg¬ 
mented  mathematical  model.  Each  segment  is  viewed  as  a  lumped 
element,  which  consists  of  six  interacting  sub-models,  and  is  con¬ 
nected  based  on  the  reactant  flow  directions  of  individual  fuel  cells. 
This  model  is  calibrated  based  on  experimental  data  of  one  single 
cell  and  the  segments  in  this  model  are  reconnected  according  to 
flow  field  patterns  of  the  other  cells.  The  distributed  properties, 
including  water  content  in  the  membrane  and  relative  humidity  in 
the  flow  channels,  are  also  investigated  by  this  model.  Modeling 
results  show  good  agreement  with  experimental  data  under  most 
working  conditions. 

We  provide  the  segmented  model  and  a  feasible  method  to  solve 
the  problem  in  which  co-flow  and  counter-flow  exist  in  the  flow 
field  design.  In  the  low  humidity  working  conditions,  current  den¬ 
sity  of  the  first  segment,  which  is  near  the  inlet,  takes  approximately 
3.9%  of  the  current  load  in  traditional  flow  field  designs.  The  value 
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can  be  improved  to  4.7%  by  placing  the  cathode  flow  direction  oppo¬ 
site  to  the  anode.  The  counter-flow  design  also  makes  water  content 
in  the  membrane  distributed  more  uniformly;  thus,  the  durability 
and  efficiency  of  the  fuel  cell  can  be  improved. 
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